A RUBIDIUM CLOCK FOR GPS 


William J. Riley 

EG&G, Inc,, Frequency and Time Department 
Salem, Massachusetts 


ABSTRACT 

The work at EGAG, Inc. on a second-source rubidium fre- 
quency standard for use in the GPS navigation satellites 
has now reached the prototype stage. This paper describes 
(1) the desip objectives and approach, (2) the more 
important design features, (3) the signal parameters and 
error budget, and (4) the early results of this program up 
to the present. 


INTRODUCTION 

Since March 1980, EG&G, Inc., has been engaged in a program to develop 
a rubidium frequency standard (RFS) for the Department of Defense, 
NAVSTAR Global Positioning System (GPS). This paper describes the RFS 
design objectives and approach as well as the early test results 
obtained up to the present prototype stage of development. 


The basic objective in the RFS development is to meet the stringent 
stability, reliability, and radiation hardness requirements of the 
Phase III NAVSTAR satellites. A summary of the general specifications 
is contained in Table 1.* The overall design approach has been to 
choose a configuration that favors the highest possible performance 
and, wherever possible, uses mature S-level quality parts. All parts 
will be subjected to extensive screening. 

The design objective of radiation hardness will be met by circuit 
design, con^ponent selection, and passive shielding. 

To ensure the highest possible performance, the design utilizes a 
discrete filter cell arrangement to permit separate control of the 
optical pumping spectrum. This allows the light shift effect to be 
independently nulled and reduces inhomogeneities in the absorption 
cell. 


♦Tables are located at the «nd of the paper. 
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The synthesizer is based on a 10.23 MHz output and uses a direct 
primary loop multiplier chain for high spectral purity. The physics 
package makes effective use of the vacuum environment for thermal 
insulation, thus reducing oven operating power. 

The early results obtained during the prototype development are 
encouraging, and the unit uses less than half the power than the 
original GPS rubidium clock (12 versus 28 watts), although it is 
somewhat larger (260 versus 170 cubic inches) and heavier (10.5 
including radiation shielding versus 9,7 pounds). The weight factor 
has only recently increased in importance, and a weight reduction 
effort is currently underway. 

It is expected that one of the prototype units will be subjected to 
extended stability testing, which will be followed by the building and 
testing of qualification units. This program should result in the 
production of the most stable rubidium clocks yet manufactured. 

SYSTEM BLOCK DIAGRAM 

A block diagram of the complete RFS system is shown in Figure 1. The 
major sections, which are discussed separately on the following pages, 
are listed below: 

(1) The physics package that serves as a frequency discriminator by 
generating an error signal whose magnitude and sense indicates 
the difference in frequency between the microwave excitation and 
the stable atomic resonance. 

(2) The supporting electronics that allow proper operation of the 
physics package. 

(3) The primary loop that produces an output signal whose frequency 
is locked to an exact submultiple of the atomic resonance. 

(4) The secondary loop that synthesizes the desired 10.23 MHz output. 

(5) The control circuits that allow remote tuning, mode selection, 
and monitoring. 

(6) The power supply. 
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Physics Package 


The physics package design is based on classical rubidium frequency 
standard principles, except for the use of a smaller TEm cavity. 
This approach gives the highest performance consistent with allowable 
size. The discrete filter cell gives zero light shift (ZLS) over a 
range of light intensity and a high S/N ratio. This permits operation 
at a relatively low light level, thus reducing temperature and rf 
power dependencies. 

The characteristics of the lamp, filter, and absorption cell are shown 
in Table 2. The lamp operates in the Kr-Rb mixed mode and is excited 
with 0.45 watt of rf power at 110 MHz. The lamp Di line emission at 
794.7 nm is shown in Figure 2a; it is free of self-reversal. Lamp 
life is aided by tight heat sinking, alkali resistant glass, high 
vacuum processing, and calorimetric measurement of rubidium fill. 

An example of calorimetric measurement is ^shown in Figure 3. A 
differential scanning calorimeter, Perkin-Elmer DSC-2C, is used to 
measure the heat energy required to melt the rubidium in a lamp. This 
allows (using the known heat of fusion of Rb) the amount of rubidium 
to be determined with a resolution of a few yG. A test program is 
underway to evaluate lamp life using this technique. 

The filter cell is operated in a separate oven whose temperature is 
adjusted for ZLS. The spectrum of the filtered light is shown in 
Figure 2b. Notice that ZLS filtration corresponds to only a slight 
attenuation of the skirt of the unwanted hyperfine optical component. 
Heavier filtration with a longer or hotter isotopic filter gives a 
stronger signal and lower filter cell TC but only at the expense of a 
large negative light shift coefficient. The filter cell buffer gas 
and pressure were chosen to give optimum signal and lowest TC under 
ZLS conditions. 

The absorption cell buffer gases are chosen for narrow linewidth and 
low TC. The mix ratio is adjusted for a slightly positive TC to 
partially compensate for the filter cell TC. The nominal fill pres- 
sure is determined by the synthesizer ratio and is sufficient to 
reduce the wall relaxation rate without excessive buffer gas colli- 
sional broadening. The absorption cell length is optimized for 
maximum signal at the chosen light intensity and temperature, the 
latter being the coolest practical value for the required upper 
ambient limit of 45®C. 
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A cross-sectional view of the physics package is shown in Figure 4. 
The three ovens are supported by a fiberglass structure with finger- 
like tabs that combine rigidity with low thermal conductance. This 
makes effective use of the vacuum environment for thermal insulation. 
Low thermal conductivity leads and low emissivity gold foil further 
reduce losses so that the entire physics package requires only 1.8 
watts of oven power under nominal operating conditions. 

It is important to avoid any condition that can distort the shape of 
the absorption line, because any such inhomogeneity will contribute to 
the rf power dependence. A lens is used to collimate the lamp output 
and provide a uniform light distribution, the filter cell is adjusted 
for ZLS, a two-section C-field coil configuration provides a uniform 
magnetic bias field, and a well bonded absorption cell with modest 
pressure shift coefficient provides low temperature gradients. All 
oven windows are sapphire for high thermal conductivity. 

Another consideration is the residual magnetic field of the oven 
heaters which, particularly for the cavity, can induce a pseudo 
temperature coefficient. Precise registration of a thin double-layer 
blanket heater and low oven supply current reduce this effect to 
2 yG/mA or 2 x 10-1V“C. 

The microwave cavity is excited with an E-probe and has slotted end 
covers that support the desired TEm mode while allowing light 
transmission without significant microwave leakage. 

The physics package includes two 0.025 inch Hipernom cylindrical 
magnetic shields. The major physics package parts are shown in 
Figure 5. 

Supporting Electronics 


Directly associated with the physics package are the temperature 
controllers, lamp exciter, photodetector preamplifier, and magnetic 
bias current supply. These circuits determine, to a large extent, the 
overall frequency stability that is achieved. 

The temperature controllers (Figure 6) are dc thermistor bridges and 
dissipative regulators, with static thermal gains of about 2000. The 
oven demand power is determined by the one hour warmup required for 
ground testing in air. Vacuum conditions not only reduce the oven 
losses but also raise the thennal gains and servo stability margins. 
Care is taken to hold thermistor self-heating constant at a fraction 
of a degree C. 
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Figure 6 — Basic temperature controller circuit 

The lamp exciter (Figure 7) is a straightforward Clapp power oscilla- 
tor with the lamp located inside the series -tuned coil. The lari^ 
network presents a range of loads as a function of lamp mode, and the 
circuit has been characterized under these various conditions to 
ensure proper operation. Starting may require the exciter to redis- 
tribute condensed rubidium inside the lamp by rf induction heating. 
Ignition takes place when sufficient voltage exists across the lamp 
coil. The running condition is stabilized against environmental 
changes and supply ripple by a current regulator circuit. 



Figure 7 — Basic lamp exciter circuit 
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The photoctetector preamplifier (Figure 8) is a two-stage configuration 
with dc and ac transimpedances of 100 kQ and 5 Wl, respectively. The 
ac gain is broadly peaked at the fundamental modulation frequency. 


O.OOIpf 



Figure 8 — Basic preamplifier circuit 


The magnetic bias supply (Figure 9) is an active current source, 
programmable under ground control with a resolution of 3 x lO’^^. 
Relay switched fixed resistors are used for the four most significant 
of 12 bits, in order to minimize tenn)erature effects on this critical 
circuit. The reference voltage for the C-field supply is derived from 
a precision zener diode located inside the oven of the secondary loop 
crystal oscillator. Provision is made for the direct use of this 
oscillator in a backup mode independent of the rubidium reference. 

ri 4 MS BITS 



Figure 9 — Basic magnetic bias supply circuit 
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Primary Loop Electronics 


The servo amplifier, crystal oscillator, and rf multiplier chain that 
complete the primary loop are shown in Figure 10. No synthesis or 
mixing is done in this loop to avoid instability caused by asymme- 
trical components on the microwave spectrum. The amount of frequency 
offset caused by a single spurious component scales directly with its 
power and inversely with its separation from the main carrier. A 
typical sensitivity is 5 x 10"13 for a component equal in power and 
5 MHz away. Direct multiplication essentially eliminates this effect. 
The most critical components are at ±18.7 kHz (from the secondary 
loop synthesizer) that are well suppressed (>-60 dBC) and symmetrical. 



Figure 10 — Primary loop block diagram 

The servo amplifier (Figure 11) consists of a main fundamental detec- 
tor/integrator channel, a second harmonic/acquisition channel, and 
the modulation reference logic. The modulation rate is derived by 
division of the 10.23 MHz output to approximately 146 Hz. The servo 
amplifier circuits, although straightforward, require care in execu- 
tion and component selection, since the system is a frequency lock 
loop and offsets directly cause frequency errors. 


INTEGRATOR 
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The crystal oscillator (made by Frequency Electronics, Inc.) uses an 
FC-type crystal which, by virtue of its excellent temperature charac- 
teristic, avoids the need for ovenization. 

The rf multiplier chain (Figure 12) consists of a phase modulator, two 
push-push doublers, and a X81 step recovery diode (SRD) multiplier. 
An ALC loop is used to maintain constant drive to the SRD nwl tipi ier. 
This approach gives high stability and spectral purity with minimum 
complexity. Particular care is taken to avoid frequency offsets 
caused by even-order modulation distortion. A relative second har- 
monic distortion level of -70 dB can cause an offset of 5 x 10"12 
and thus must be constant to a few percent to obtain good frequency 
stability. A pure modulation waveform is generated by passive inte- 
gration of a precision squarewave, and a highly linear phase modulator 
characteristic is obtained by applying small excursions to a hyper- 
abrupt varactor diode biased at an "S" shaped inflection point. A 
more subtle factor is the distortion caused by AM and subsequent 
AM-to-PM conversion in the multiplier chain. A simple RC phase 
shifter circuit introduces equal amounts of AM and PM and thus can 
cause considerable distortion. The all -pass phase modulator configur- 
ation shown in Figure 12 provides about 26 dB of AM rejection and thus 
greatly reduces this problem. 
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Figure 12 — Basic rf multiplier chain and circuit 
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The SRD multiplier (made by Frequency Sources, Inc.) uses a low 
capacitance diode in a shunt mode with external dc bias. A tetn)era- 
ture stabilily of 0.01 dB/^C is typical for the **50 uW microwave 
output which can be adjusted to the optimum level by an internal PIN 
diode attenuator. 

Secondary Loop and Synthesizer 

The 10.23 MHz output is derived from a precision ovenized crystal 
oscillator (made by Frequency Electronics, Inc.) that is phase locked 
to the primary loop by the synthesizer. (See Figure 13.) The +18 dBm 
output is obtained via an output amplifier that has sufficient tank 
circuit Q to maintain a phase-continuous output under transient 
radiation. Sufficient isolation is provided by the divider input 
circuits to keep the output spectrum free (-100 dBC) of synthesizer 
spuri i . 


w 18.7 KHz 



AMP. 


Figure 13 - Synthesizer block diagram 

The synthesizer consists of a pair of dividers that divide the primary 
and secondary loop frequencies down to their common factor of **18.7 
kHz. A phase detector and loop filter generates a control voltage for 
the 10.23 MHz crystal oscillator. Rapid lock acquisition is accom- 
plished by logic that presets the dividers to the proper phase rela- 
tionship. This also provides fast recovery from a radiation-induced 
upset. 

Power Supply 

The power supply consists of an input filter, a dc/dc converter, and 
two linear regulators as shown in Figure 1. The dc/dc converter 
(Figure 14) uses a single-ended flyback configuration that provides 
both dc isolation and regulation. Its +28 V and +5 V outputs are used 
directly as heater and logic supplies respectively while the +16 V 
outputs are further stabilized by precision linear regulators before 
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supplying the analog circuits. Under normal operating conditions In 
vacuum, the RFS circuits consume under 10 watts. The total dc Input 
Is expected to be 28 watts demand during warmup and 12 watts steady 
state. 



Figure 14 — Basic power supply circuit 


Packaging 

The overall RFS package, shown In Figure 15, consists of an aluminum 
structure with aluminum and Hipernom covers for magnetic, EMI, and 
radiation shielding. The hollow center of the structure holds the 
physics package, secondary crystal oscillator, the lamp exciter, 
preamplifier, and SRD multiplier. The other electronic circuits are 
on 13 circuit boards located Inside compartments that comprise the 
outside walls of the structure. Additional local radiation shielding 
Is used as required. This approach meets all thermal, electrical, 
mechanical, magnetic, and radiation requirements as well as allowing 
good access to most circuit boards. 



Figure 15 — RFS outline drawing 
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Signal Parameters 


The operating conditions and signal parameters of the physics package 
are summarized in Table 3. In general, these values indicate a good 
signal and exceptionally low sensitivities. 

The S/N ratio of a passive rubidium frequency standard is limited by 
the shot noise of the dc photodetector current, Iq, which sets a 
level of white frequency noise and gives a time-domain frequency 
stability that varies as where t is the averaging time. 

The frequency stability can be estimated simply fr(xn the values of 
Iq and the discriminator signal , 0. For Iq * 82 yA, the shot noise 
current, in. is: 


ip = 5.1 pA/-0iz 


This establishes the white freouency noise leveVin relation to the D 
value of 240 pA per 1 x 10“10, giving a fractional frequency noise 
power spectral density of: 

Sy(f) = ho = {— X 10-10)^ = 4.5 X 10-24(Hz)-i 


This results in a time-domain stability of: 

Oyh) = 1.5 X 10-12(t)-1/2 

which may be converted to a frequency-domain value of: 


^(f) = COy{x)]2 Vq^X f"2 


where Vq is the carrier frequency of 10.23 MHz. Then, 


i^(f) = 2.4 X 10-10 f-2 and;C(l) = -96 dBC/Hz 


The measured noise level is 5.5 pA-^liz rather than the 5.1 pA/yiiz 
value expected from shot noise only. The noise level is 0.7 pA/^Hz 
without light, also indicating a negligible contribution from the 
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photodetector and pre-amplifier. This is possible primarily because 
of the use of a very high quality EG&G custom photodetector that has a 
high leakage resistance (>50KO ) at the 65®C operating temperature. 

The actual noise results in an estimated time-domain stability of 1.6 
X 10"12 x-1/2. The observed value is 4 x i-l/2 uncorrected for 

the (unknown) contribution of the reference. 

The S/N ratio can be defined as the rms value of the maximum fundamen- 
tal error signal current divided by the rms noise current in a 1 Hz 
bandwidth. 

Values of over 90 dB have been obtained with this design by operating 
the lamp in the more intense red (all Rb) mode, and over 80 dB in the 
mixed (Kr-Rb) mode. However, higher light intensity also results in 
larger temperature coefficients, rf power shift, and other sensitivi- 
ties. For the GPS application, it is considered a better tradeoff to 
favor long-term and environmental stability by operating at modest 
light intensity where the S/N ratio is about 75 dB. 

ERROR BUDGET 

The RFS frequency is influenced by a number of factors that can affect 
its stability. Some of these factors involve dependencies of the 
physics package and others are a result of electronic circuit sensi- 
tivities. The largest of these factors are listed in Table 4. Also 
shown in Table 4 are the resulting individual temperature coefficients 
and the total worst-case and rms values expected. These factors also 
govern the aging rate and radiation sensitivity of the unit. 

TEST RESULTS AND CONCLUSIONS 

Tests of a breadboard version of the RFS have given encouraging 
results. The results of stability tests conducted on the primary loop 
with the physics package in a thermovac chamber and most electronic 
circuits on the bench in air are shown in Figure 16. The measured 
time-domain stability is within the goal specification. The results 
are uncorrected for the contribution of the reference and are probably 
limited by the reference stability in the medium term region. The 
region near x = 1 second will be determined by the secondaiw loop 
crystal oscillator, which has a measured stability of 4 x The 

RFS stabili^ in the medium term region is at least as low as 
4 X 10"12 x-i/2 and reaches below 1 x 10"13 between x= lO^ and lO^ 
seconds. The Rb reference is better in the medium term region, but 
the Cs reference has been found to be better at 10^ seconds and 
longer, probably because of the barometric sensitivity of the rubidium 
reference unit. The GPS RFS is immune to this effect in the thermovac 
chamber. A drift rate of under 1 x 10"13 day and a physics 
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package temperature coefficient of 1 x has been obtained. 

It is reasonable to assume that the 10® second stability will be 
even better in the prototype version, when all circuits are packaged 
and operated within the controlled thermovac environment. Overall, 
the RFS program is on schedule and all objectives are being met. 



Figure 16 — Time-domain stability of the breadboard 
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Table 1 — Condensation of specifications for the Phase III GPS RFS 


Output 

f =10.23 MHz nominal, +18 dBm nominal 

Operating Temperature 

+ 20 to +45*C 

Frequency Stability 
{Figure 16) 

®y(T) = 1 X 10“11 for X = 1 and 10 
seconds decreasing as^ to goal of 
lx 10“ 13 at 105 seconds 

Phase Noise 

X(l) = -83 dBC/Hz maximum, decreasing 
at 10 dB/decade to 5 kHz 

Spurious 

-100 dBC maximum, 5 kHz to 20 MHz from 
carrier 

Temperature Coefficient 

i 1 X 10-12/°c 

Drift 

■1 

CO 

1 

o 

>< 

VI 

<1-1 

j 

Voltage Sensitivity 

^ < 1 X 10-12/volt 
f 

Mag Susceptibility 

^ < 1 X 10-12/3 Gauss 
f 

Warmup Time 

1 hour 

Reliability 

73 % for 7.5 years' operation 
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Table 2 — Lamp and cell characteristics 



Lamp 

Filter Cell 

Absorption Cell 

Glass (Corning #) 

1720 

7070 

7070 

Outside Diameter* 
inches 

0.31 

1.09 

1.09 

Inside Length, 
i nches 

0.54 

0.18 

0.45 

Isotope 

87 

85 

87 

Buffer Gas 

Kr 

Ar 

= 0.731 

Pressure, Torr 

7 

...... 

120 

14.2* 

Operating 
Temperature, ®C 

120 

85 

65 

Temperature Coefficient 
pplOll/“C 

~o 

-6 

+3 


♦Offset = +2.78 kHz. 
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Table 3 — Breadboard RFS operating characteristics 


Si gnal 


Second Harmonic Signal: 

45 mV rras* 

Maximum Fundamental Signal: 

165 mV rms* 

Light: 

8.23 Vdc* 

Discriminator Slope: 

1.2 mV per 1 x IQ-IO* 

Linewidth: 

252 Hzt 

Sensitivities 


Light Shift: 

~0 + 1 X 10“12/% 

Lamp Oven TC: 

~0 ± 3 X 10-12/oe 

Filter Cell TC: 

6 X io-n/®c 

Absorption Cell TC: 

3 X lO-U/^C 

RF Power Shift 

ss± 2 X 10-12/dB 

Lamp Exciter TC: 

-0 ± 5 X 10-1 V*C 

Noise 


Shot Noise: 

5.1 pA/^llz 

Total Noise: 

5.5 pA/'^/Hz 

Si gnal : 

33 nA rms 

(S/N): 

76 dB in 1 Hz BW 

Line Q: 

16 X 106 

Predi cted Stabi 1 i ty : 

ay{T) = 1.6 X 10-12 t-1/2 


*At preanq) output with Zj(ac) = 5 MS, Zj(dc) = 100 Kfi. 

"^Full width between inflection points with normal light and rf. 
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Tab! e 4 — Frequency error budget 


Factor 

Sensi ti vi ty 

TC, pplOlV^C 

Magnetic Bias Current 

+ 3 X 10"lVppm 

± 6.0 

Lamp Exciter TC 

± 5 X 10-1V*’C 

± 4,0 

Filter Cell TC 

- 6 X lO-lV’C 

- 3,0 

d^/dt Effects 


± 2.5 

RF Level Shift 

+ 2 X 10-1 2/dB 

±2.0 

Residual Heater Field 

± 5 X 10-lVmA 

±2.0 

Absorption Cell TC 

+ 3 X 10-11/"C 

+ 1.5 

VCXO Dynamic TC 

— 

± 1.5 

Servo Integrator 

i 1 X 10-14/“C 

± 1.0 

Power Supply 

t 1 X 10-12/% 

± 1.0 

Modulation Deviation 

± 1 X 10-13/% 

± 1.0 

Miscellaneous 

— - 

± 1.0 

Total 

Worst-case 23,5 



rms 9.2 


629 









